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One of the major open issues in neutrino physics is the possible existence of CP violation in 
the neutrino sector. Such an observation would have an important impact in various domains of 
physics, from high energy physics to cosmology. Its search requires future accelerator neutrino 
facilities producing intense and pure neutrino beams such as "beta-beams". Here we review the 
different beta-beam scenarios proposed so far and discuss the present status, with a particular 
emphasis on the original baseline scenario and its feasibility. Alternative strategies for the CP 
violation search are to be pursued as well. A possibility is to search for CP violation effects in 
astrophysical environments. Here we present recent analytical and numerical results obtained in 
the context of core-collapse supernovae. In particular, we point out the conditions under which 
there can be CP violating effects in dense media and show numerical results on the supernova 
(anti-)neutrino fluxes and on the electron fraction, relevant for the r-process nucleosynthesis. 
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1 Introduction 

On 

O ■ In the last decade a series of experiments have provided us with essential information on neutrino 
properties, giving clear evidence for physics beyond the Standard Model. Neutrino oscillations first 
^ . proposed by B. Pontecorvo [lj, are a firmly established phenomenon and most of the oscillation param- 
eters have now been measured with accuracy. The ensemble of the present experimental results can be 
interpreted in the context of three active flavours. In particular, the characteristic L/E u dependence^ 
has been identified as well both by the Super-Kamiokande and the Kamland experiments. The oscil- 
lation phenomenon implies that the flavour and the mass basis are related by a unitary matrix, the 
Maki-Nakagawa-Sakata-Pontecorvo matrix [2] : 

fl 

e ia 
\0 e i/3 

with Cij =cos9ij (sij =sm9ij) and #i 2 , #23 and #13 the three neutrino mixing angles. Two angles and 
two squared mass differences have been precisely measured, namely 5m\ 2 = 8 x lCT 5 eV 2 , |#mf 3 | = 
3 x 10 _3 eV 2 , sin 2 2# 12 = 0.83 and sin 2 2# 23 = 1 thanks to the solar and reactor experiments on one hand 
and to the athmospheric ones on the other [3]. 




1 L is the distance between the neutrino source and the detector, E v is the neutrino energy. 
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Four important neutrino properties remain unknown : the value of the third mixing angle, the ab- 
solute neutrino mass scale and hierarchy, the (Dirac versus Majorana) neutrino nature and the possible 
existence of CP violation in the lepton sector. Upcoming experiments and projects that are at present 
at the R&D level will address these questions. For the value of #13 we only have an upper limit given 
by the Chooz experiment [I]. In the near future three reactor experiments - Double-Chooz [5], RENO 
[6J, Daya Bay [7] - and the first super-beams experiments (T2K and NOi/A) will be able to measure its 
value if sin 2 2#i3 < 0.02 |8J. Note that the combination of the available experimental data gives indica- 
tion that #13 might be close to the present Chooz limit [9j [TDj . The second unknown parameter, relevant 
for oscillation experiments, is the Dirac phase 5. Only if the third neutrino mixing angle is different 
from zero there can be CP violating effects coming from the Dirac phase. Indeed for a non-zero 5 the 
MNSP matrix becomes complex, introducing a difference between neutrinos and anti-neutrinos. The 
matrix in Eq.([T]) is completely specified if the Majorana phases a,/3 are also determined through the 
measurement of neutrinoless double-beta decay that addresses the crucial issue of the neutrino nature 
and of the (effective electron) neutrino mass. It is clear that the upcoming experiments will tell us if 
#13 is close to the present Chooz limit or if it is (very) small. This will set our future strategy for CP 
searches with accelerator experiments. 

As is well known the reactor experiments are not sensitive to the Dirac phase since they measure 
electron (anti-)neutrino disappearance and the electron neutrino survival probability in vacuum does not 
depend on 5. To this aim one needs appearance experiments such as super-beam ones where oscillations 
of to v e are looked for. As pointed out in [S], if sin 2 2# 13 < 0.02 the upgrades of the first generation 
super-beams (NOz/A and T2K) have a CP discovery reach for a significant fraction of 5 values at 90% 
CL while a 3<r discovery can be found only for a tiny fraction of the phase values. Therefore, unless we 
are in a very lucky situation in which the S value is fixed at one of the few values that can be covered 
with these facilities, the measurement of the Dirac phase will require intense neutrino beams, based 
on new concepts, - the neutrino factory or the "beta-beams" -, or the second generation super-beams. 
The neutrino factory exploits neutrino beams produced by the decay of stored muons and sent to far 
detectors. The beta-beam concept, first proposed by P. Zucchelli [IT] , uses the beta-decay of boosted 
radioactive ions. The goal of these facilities is to search for small #13 values, determine the Dirac phase 
and identify the mass hierarchy. Here we discuss the different beta-beams scenarios proposed - the 
original, the high energy scenarios [12] as well as the electron capture variant [32] - and discuss some 
of the feasibility aspects. (For a review of all scenarios see [II].) We will also mention low energy 
beta-beams [15J, with neutrino energies in the 100 MeV energy range. The availability of such beams 
would broaden the physics case since they offer the unique opportunity to perform neutrino interaction 
studies of interest for nuclear structure, for the study of fundamental interactions and for the physics 
of core-collapse supernovae. 

In paraller with the feasibility and physics reach studies of long term facilities such as neutrino 
factories and beta-beams, it is important to explore complementary avenues e.g. by searching for 
indirect CP violation effects in astrophysical environments. So far few attempts have been performed in 
this direction. In [16] the authors have studied the effects of the Dirac phase in a star like our Sun. In 
particular, it has been shown that the electron neutrino survival probability does not depend on the CP 
phase even in matter, while the effects coming from radiative corrections are expected to be extremely 
small. In [T7J CP effects in core-collapse supernovae are considered in the extreme density limit. A 
somewhat different perspective is followed in [18] where it is pointed out that the combination of an 
early CP measurement with Earth based experiments and with neutrinos from astrophysical sources 
can significantly improve the CP discovery reach. Recently we have been exploring the conditions under 
which there can be CP violating effects in supernovae [12]. Exact analytical results have been obtained, 
for the first time, that are valid for any density profile. Besides, extensive numerical calculations have 
been performed quantifying possible effects on the neutrino fluxes in the star, in an observatory on 
Earth and on the r-process nucleosynthesis. These results have been extended in a following work [20] 
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Figure 1: The beta-beam scenario based at CERN using existing accelerators (the PS and 
the SPS), in conjunction with the EURISOL facility, as first proposed by Zucchelli The 
low energy accelerator infrastructures are according the feasibility study performed within 
the EURISOL (European Isotope Separator On-Line Radioactive Ion Beam Facility) Design 
Study (6th Framework Programme) [21]. The neutrino beams produced in the storage ring 
straight sections are sent to a large-size detector, such as MEMPHYS, located 130 km from 
CERN. 



including in particular the neutrino-neutrino interaction. The works in [T9l [20] have set the basis for 
the exploration of CP violating effects in dense media. 

2 Beta-beams 

The CP violation search will require neutrino beams with the highest possible intensities, tiny intrinsic 
backgrounds and a very good control of systematic errors. The beta-beam concept has three important 
advantages. The neutrino beams are pure in flavour since only electron neutrinos or anti-neutrinos can 
be produced, depending on the ion that decays through (3 + or (3~ . This means that there is no beam 
related background. The neutrino intensity and energy spectrum is precisely known, since the number 
of ions is perfectly controlled. The flux emittance is inversely proportional to the Lorentz boost factor 
7- 

The beta-beam scenarios proposed so far can be summarized as follows [14J : 

• the standard beta-beam scenario : In the original baseline scenario [11], the beta-beam 
facility is hosted at CERN (Figured]). The ions are produced through the ISOL technique, bunched 
and accelerated first at several hundred MeV per nucleon, then injected in the PS and SPS to 
reach 7 = 100 (Figure [2]). Finally the ions are stored in a storage ring with long straight sections 
(2.5 km for a total length of 7 km) that point to a (large size) water Cerenkov detector, located 
in an (upgraded) Underground Frejus laboratory. The search for CP effects is performed through 
the comparison of v e — > z/ M versus v e — » oscillations. Note that if one sends a super-beam to the 
same detector, T breaking can also be studied by measuring —>■ v e oscillations as well as CPT 
breaking through the comparison of v e — > and —> v e oscillations. The discovery reach for CP 
violation of the original scenario has been first investigated in [24J and in detail in [23] . With such 
a baseline one has a sensitivity on the third neutrino mixing angle of sin 2 #i3 > 5 • 10 -4 (3<r for 75 % 
of the values of 5) and has a CP violation discovery reach at 3a for sin 2 #i3 values down to 7- 10~ 4 
(Figure [3]). While the baseline is too short to exploit the matter effects in order to determine the 
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Figure 2: Left : Comparison of neutrino fluxes from a super-beam (SPL) and the standard 
beta-beam with 7 = 100 [22]. Right : Electron anti-neutrino fluxes from low energy beta- 
beams, produced via the decay of 6 He ions boosted to 7 = 7 (dotted line), 7 = 10 (broken 
line) and 7 = 12 (dash-dotted line) [23] . 



mass hierarchy, this goal can be achieved by using a combination of the neutrino beams with the 
athmospheric ones as first proposed in |26j and used in [25] to show that the octant and the mass 
hierarchy identification can be determined with the standard beta-beam (Figure HJ). 

• the low-energy beta-beam [14] : This corresponds to ions running at 7 = 5 — 15 to pro- 
duce neutrinos in the 100 MeV energy range (Figure [2]). Such beams would allow : i) neutrino 
scattering measurements on nuclei, to improve our knowledge of the isospin and spin-isospin nu- 
clear response [HI [271 ESI EH] for various timely applications (for example to costrain neutrinoless 
double-beta decay calculations [30J, r-process calculations and the detection of (relic) supernova 
neutrinos; ii) searches of physics beyond the Standard Model (through e.g. the measurement 
of the Weinberg angle [23J, a new Conserved Vector Current hypothesis test [31 J, of coherent 
neutrino- nucleus scattering measurement [32], searches of non-standard contributions [33]); iii) 
core-collapse supernova physics studies [HI EH [35] • This option might require a devoted storage 
ring [27] or could use the low energy neutrino beam component of the standard beta-beam by 
putting one/two detectors at off- axis [36] . 

• the high energy beta-beam [12] : In the first proposal the ions run at higher 7 (about 350) 
and the baseline is increased accordingly to about 700 km. Such a long baseline is sensitive to the 
neutrino hierarchy through matter effects. Note that increasing the boost factor require a bigger 
storage ring since the length increases linearly with 7. On the other hand a higher 7 produces 
higher rates in the detector for the same ion intensity. The work in [12] has stimulated numerous 
studies (see e.g. [37J [381 EH STJ1 HU H2])- For example in [13] the magic baseline from CERN 
to the India-based Neutrino Observatory (7152 km) is investigated. At such a distance the v e 
to Up survival probability has no dependence on 6 and allows to measure the neutrino hierarchy 
without any degenerate solution. A beta-beam experiment from FNAL to DUSEL is proposed 
in [H]. Boosting the v e and v e to the maximum energy possible with the Tevatron one could 
exploit an interesting combination of the first and the second maximum that helps reducing the 
parameter degeneracies. 

• the electron capture option [13] : Here the ions decay by electron capture instead of beta- 
decay, producing electron neutrinos only. The advantage of such beams is that the neutrinos 
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Figure 3: Sensitivity on 5 at 3<r for the (7 = 100) beta-beam (/3B), a super-beam (SPL), 
and T2HK (see text) as a function of the true value of sin 2 2#i3. The width of the bands 
corresponds to values with 2% to 5% systematical errors. The dashed curves correspond 
to the combination of the beta-beam and the super-beam with 10 yrs measurement time 
each and 2% systematical error. For parameter values inside the ellipse-shaped curves CP 
conserving values of 5cp can be excluded at 3<r (A^ 2 > 9) [25] . 



are monoenergetic. On the other hand this option requires acceleration and storage of partially 
stripped ions, which is quite challenging from the technical point of view. Candidate ions have 
rather long lifetimes like for example 150 Dy. Promising sensitivities on CP can be obtained by 
running the ions at two different 7 tuned to the first and second oscillation maximum [15]. A 
combination of a beta- and electron-capture beam is investigated in [16] considering different 
baselines. 

The physics reach of the various beta-beam scenarios has been investigated in great detail in the 
last years. The studies performed can be classified in two categories, based on what one can call 
the "conservative" or "optimized" attitude. In the former the ion intensities and boosts are obtained 
extrapolating well known technologies and existing accelerator infrastructures (the present PS and SPS 
at CERN for example). In the latter the ion intensities and boosts are treated as "free" parameters, 
with the aim of exploring the conditions to achieve optimal sensitivities. Figure [5] presents a comparison 
on the CP discovery reach of the future long-term facilities [47] . For each option the figure shows the 
sensitivity corresponding both to the conservative and to the optimized options. The T2HK (SPL) 
cases correspond to a 4 MW proton driver, with the neutrino beams firing to a megaton class (440 kt) 
water Cerenkov detector located at 295 (130) km distance. The running time is 2 (8) years for neutrinos 
(-neutrinos). The conservative option for the beta-beam is the CERN-Frejus scenario with 7 = 100, 
while the optimized one is with a 7 = 350 and a baseline of 700 km. In both cases a 440 kt water 
Cerenkov detector is considered and the running time is 5 (5) years with 2.9 10 18 (1.1 10 18 ) 6 He ( 18 Ne) 
per year. For the neutrino factory, the conservative (optimised) setup uses 10 21 muon decays per year, 
with a stored muon-beam energy of 50 (20) GeV, The beams point to a (two) 50 kton detectors located 
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Figure 4: The mass hierarchy and the octant degeneracy can be identified in the standard 
beta-beam scenario by exploiting the synergy with athmospheric neutrinos. Three cases are 
shown : the 7 = 100 beta-beam, the super-beam (SPL) at CERN and T2K phase II to 
Hyper-K. The figure shows the allowed regions in sin 2 2# 13 and #cp for accelerator data alone 
(contour lines) and accelerator plus atmospheric data combined (colored regions). Solution 
for the true/wrong mass hierarchy (octant of 6*23) are indicated with H tr / wr (O tr/wr ). The true 
parameter values are <5cp — — 0.857T, sin 2 26> 13 = 0.03, sin 2 ^ = 0.6 [25] . 



at 4000 (4000 and 7500) km. The running time is 4 (5) years with fi + and fi~ . As it can be seen from 
Figure [5] if sin 2 2# 13 >0.02 the three options have very similar sensitivities on 5. For values of 5 10~ 4 < 
sin 2 26' 13 < 10~ 2 the super-beams are outperformed by the beta-beam and the neutrino factory. Only 
the optimised neutrino factory can reach values smaller than sin 2 2# 13 < 5 10 -4 [47], while the optimised 
beta-beam and the conservative neutrino factory options have a comparable performance. 

A first feasibility study of the original beta-beam scenario is performed in [48], while a detailed 
investigation has just been completed within the EURISOL Design Study [_T] . An important step 
forward has been made. Indeed most aspects appear under control (e.g. bunching of the ions, the 
stacking method in the storage and contamination issues). Reaching the required ion intensities will 
need further investigations [19]. Two sets of ions have been discussed so far, namely 6 He and 18 Ne 
on one hand and 8 Li and 8 B on the other. The production of 6 He and 18 Ne through the standard 
ISOL technique as well as the direct production method have been studied within the EURISOL DS. 
The production of 8 Li and 8 B with a storage ring method has been recently proposed [50J and will be 
investigated within the EUROnu initiative. 

3 CP violation and astrophysics 

3.1 The conditions for CP violation in core-collapse supernovae 

Our understanding of neutrino propagation in massive stars has been revolutionized in the last few 
years after the use of temporally evolving density profiles with shock waves and the inclusion of the 
neutrino-neutrino interaction. This is very different from the solar case where the neutrino propagation 
is well understood in terms of the standard MSW paradigm. Indeed new phenomena have emerged, such 
as multiple resonances and phases effects on one hand, and of collective effects on the other. Besides, 
the study of the impact of turbulence is just at its beginning. (For a review see |51j). 

Recently we have investigated possible CP violating effects coming from the Dirac phase in super- 
novae. Here we describe the conditions under which there can be CP effects. We use the procedure 
established in [19] and sketch the main steps of the demonstrations : 



6 



Sensitivity to LCPV 




10~ 5 10 " 10"' 10"' 10 1 

True value of sin 2 2613 



Figure 5: The CP discovery reach for various proposed long-term accelerator facilities : 
super-beams (T2HK in yellow, SPL in orange, WBB in green), beta-beams (light green 
band) and a neutrino factory (blue band). The sensitivities are shown for the fraction of 
all possible values of the true value of the CP phase 5. Each band is determined by two 
different options for each facility, a conservative one giving the lower end while an optimised 
one gives the upper end (see text) [27] . 



within the standard MSW effect at tree level 



with the neutrino-neutrino interaction term with/without the V^ T potential coming from the 
radiative contributions 



3.1.1 The neutrino evolution equations 

Let us consider the three flavour neutrino evolution equations in matter : 
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where V c (x) = y/2GpN e (x) for the charged and V n {x) = — ^GpN n (x) for the neutral currents. Since 
V n only contributes an overall phase to the neutrino evolution it can be ignored. The V^ T term is due to 
radiative corrections and is such that V^ T /V C ~ 10~ 5 [52]. The last term of the Hamiltonian corresponds 
to the neutrino-neutrino contribution and is given by : 



h w = V2G F EE/ Mq')(i - q ■ %) dn c 

1/ r, J 



(3) 



a v a ,v. 



where Gp is the Fermi coupling constant, p = p Ua (—p* Ua ) is the density matrix for neutrinos (anti- 
neutrinos), q (q') denotes the momentum of the neutrino of interest (background neutrino) and dn a is 
the differential number density. In the single-angle approximation, that assumes that the neutrinos are 
all emitted with the same angle, i.e. p(q) = p(q), Eq.Q reduces to 



V2Gp 
lixRl 



(4) 



with the geometrical factor D(r/R u ), where the radius of the neutrino sphere is R u , and L Ua (q) are the 
neutrino fluxes at the neutrinosphere, which can be taken Fermi-Dirac or power-law distributions. 
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3.1.2 The conditions for CP violating effects 
The standard MSW case 

To understand the origin of the CP effects, it is useful to rewrite Eq.(l2"l) in the T 2 3 rotated basis : 

*„ = cos023*/*-sin0 23 * T , (5) 

* T = Sin 023*^ + 008 023*^ (6) 

If one neglects the neutrino-neutrino interacation, after some calculations one can show that the whole 
5 dependence of the Hamiltonian can be factorized out as 

U{6) = S j U{5 = 0)5 <-> H{5) = S j H{5 = 0)S (7) 

where the CP dependence is in the unitary matrix 
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From Eq.([7|) one can show that the following exact relations on the probabilities hold 

P{y e ^v e ,5^0) = P{u e v„ 6 = 0). (9) 

and 

-> u e , 5 0) + P(u T ^v e ,5^0) = ^u e ,5 = 0) + P(u T ^u e ,5 = 0). (10) 

The first relation was first shown in [16] with a different procedure, in the context of solar neutrinos. 
Eq. fllOl) is new and has important implications on the neutrino fluxes in a core-collapse supernova. We 
emphasize that our relations Eqs. (I91 TT01) are exact and valid for any density profile. In a following work 
[53] . these and new relations among the probabilities have been obtained using the adiabatic basis. 
The consequence on the electron (anti-) neutrino fluxes can be immediately seen from 

<f>„.{6) = L Ue P{v e -> v e ) + L Ux {P(^ -> v e ) + P(v T -> v e )) (11) 

where v x = or u T , which holds if and only if if the and v T neutrino fluxes are exactly equal at the 
neutrinosphere. In this case <fi Ue (S) and </>p e (<5) do not depend on 5. 
The case with the neutrino-neutrino interaction 

To show the origin of CP effects in this case, it is better to start from the Liouville-Von Neumann 
equation for the density matrix (% = 1) : 

= [UH vac U^ + H m + H w {5\ Pv JS)\, (12) 

with H vac and H m the usual vacuum and matter contributions. To prove our result (that the CP- 
violating phase can be factorized out of the total Hamiltonian which includes H uu ), one has to rotate 
in the T23 basis : 

. dSp^)S_ = ^^^^^t + Hm + sh vv (5)S\ Sp^(5)S% (13) 



We now show that Eq.Q is indeed satisfied for the total Hamiltonian in Eq. (jl"2]) including the non- 
linear H uu term of Eq.([3]). Let us now consider the evolution equation of the linear combination 
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Figure 6: Ratios of the electron fraction for 5 = 180° compared to 5 = 0°, as a function 
of the distance from the neutron-star surface. The initial u^, v T fluxes have temperatures 
which differs by 1 MeV (see text). The results correspond to the normal hierarchy and 
sin 2 2# 13 = 0.19 [19]. 



L Ua Sp Ua {cL 1 5)S' { at a given momentum q. At the initial time, this quantity reads, in the T 23 basis 

as : 



One immediately sees that this quantity does not depend on S if and only if L v = L Vt . Moreover, 

one can show that since the total Hamiltonian of Eq. (1131) is independent of 5 at initial time (for H vu 
under the condition L v = L Ut ), this is true at any time, by recurrence from the Liouville-Von Neumann 

equation Eq. ffl2l) . The same applies to the anti-neutrinos. One can see that, since H uu is proportional 
to the quantity in Eq. f[T4"j) (both for neutrinos and for anti-neutrinos) this implies that the factorization 



Hence Eqs.(l9l- (jT0l) are valid for the total Hamiltonian with the vv term as well. Note that this result is 
true independently on the angular assumption made for the neutrino-neutrino contribution (single-angle 
versus multi-angle). 

General conditions for CP violating effects 

From the above argument it is clear that there are at least two possible sources of CP effects in dense 
media : 1) any contribution to the Hamiltonian that breaks the factorization Eq.([7|) in which case 
Eqs.fl MTOl do not hold any more; 2) any physical effect that makes the [i and r neutrino fluxes at the 
neutrinosphere differ, in which case the relation Eq. fllOl) does not hold any more. It is easy to show that 
conditions 1) is for example fulfilled when (standard or non-standard) radiative corrections are included. 
As a consequence the electron neutrino survival probability becomes dependent on 5, breaking Eq.(j9]). 
Such effect that is in principle very small, can be amplified in presence of the non-linear vv contribution. 
Physics beyond the Standard Model, such as Flavour Changing Neutral Currents, can also break the 
condition 1) and 2) at the same time. This can induce indirect CP effects in supernovae. 

3.1.3 CP effects in core-collapse supernovae : Numerical results 

We have performed numerous calculations of the CP effects on the neutrino fluxes within a supernova 
[T9l [20] . The input parameters used are the oscillation best fit values and supernova density profiles 




(14) 



in Eq.flT]) holds 



H w {8) = SH VV {5 = 0)&. 



(15) 
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Figure 7: Ratios of the v e fluxes for a CP violating phase 5 = 180° over 5 = 0° as a function 
of neutrino energy, at 1000 km within the star, without the vv interaction and radiative 
corrections (dotted), with the vv interaction (dashed), with the vv interaction and radiative 
corrections taking different /i and r fluxes (solid). The case where the vv interaction and 
radiative corrections are taken into account with v^ and v T fluxes equal (dot-dashed) is also 
shown. The results correspond to an inverted hierarchy and a small 6*i 3 [2"U] . 



that fit supernova simulations. We find significant CP effects on the muon and tau neutrino fluxes. 
However such effects do not have any impact since they disappear when one sums up the electron, muon 
and tau neutrino flux contributions since they all interact through neutral current. 

Concerning the possible CP effects on the nucleosynthesis of heavy elements (the r-process), Figure 
[6] shows the effects on the electron fraction that governs the neutron to proton ratio, a key parameter 
for the r-process abundance calculations. One sees that for this observable the effects appear to be very 
small, although the effect on the electron (anti-)neutrino fluxes are of the order of several percent. The 
CP effect in an observatory on Earth comes out to be too small to be detectable. Such effects appear 
at high neutrino energies in electron anti-neutrino scattering on protons, the main detection channel of 
water Cerenkov and scintillator detectors. Note that the size of the modifications due to the 5 phase 
does not depend on the size of the observatory [TjJ] . These results have been obtained by including the 
coupling of neutrinos with matter only (standard MSW). Further calculations should be performed to 
see how the CP effect on the electron fraction as well as in an observatory on Earth can be modified 
when one includes the neutrino-neutrino interaction. 

Figure [7J shows the CP effects on the neutrino fluxes in the star when the radiative corrections 
as well as the neutrino-neutrino interaction contribution are included as well [211] ■ Several differences 
appear when this term is included. In particular, the synchronization regime that is present in the inner 
part of the star tends to freeze the CP effects that only appear during the bipolar oscillations regime. 
Contrarily to the calculations without the vv term, it is the low energy neutrinos that are affected by 
the CP phase instead of the high energy ones. The neutrino fluxes in the ten MeV energy range are 
modified by approximately 5 — 10% when the v^ and the v T neutrino luminosities at the neutrinosphere 
differ by a few percent (as an example). Such effects are also present if the mu and tau neutrino fluxes 
at the neutrinosphere are taken equal if radiative corrections are taken into account. Indeed the size of 
the CP effects is usually at the level of ~ 10~ 4 — 10" 5 when the radiative corrections only are included. 
This is in agreement with what was found in [T6] in the context of solar neutrinos. It is the non-linearity 
of the neutrino evolution equations that amplifies such effect at the level of a few percent, as shown in 
Figure [7J [20] . Note that non-standard contributions (such as those from SUSY models) could make the 
radiative corrections much larger [51]. How much the CP effects can be amplified in this case requires 
further investigation. 



10 



4 Conclusion 



Crucial open questions remain in neutrino physics, among which the value of the third neutrino mixing 
angle, the neutrino nature, the mass scale and the hierarrchy, and the possible existence of CP violation 
in the lepton sector. If sin 2 2#i3 > 0.02 the third mixing angle will be measured by next reactor and the 
first super-beam experiments. However, unless nature has fixed the CP phase to one of the few values 
that can be covered with the upgrades of such super-beam facilities, the search for the CP violating 
Dirac phase will require long-term accelerator facilities, such as the second generation super-beams, the 
neutrino factory and the beta-beams. As far as beta-beams are concerned, many baseline scenarios 
have been proposed so far, for which the physics reach is under intense investigation. The CERN-Frejus 
baseline has now been studied in depth both from the point of view of the physics reach and of the 
feasibility, thanks to the EURISOL Design Study just ended (FP6). Further investigation of the ion 
intensity with the ISOL technique, the direct production and the storage ring methods need to be done/ 
The standard baseline is very attractive if #13 is large or not too small. In particular, it can pin down 
if sin 2 #i3 > 5 • 10 -4 (3a for 75 % of the values of 5) and has a CP violation discovery reach at 3a for 
sin 2 # 13 values down to 7 • 10~ 4 . The octant degeneracy and the mass hierarchy can be identified as well 
by exploiting the synergy with athmospheric neutrinos in the same large size detector. 

It is important to pursue alternative searches for CP violation through the exploration of indirect 
effects e.g. in astrophysical contexts. Recently it has been demonstrated that there can be CP violating 
effects in core-collapse supernovae. These can emerge for example from radiative corrections, from 
physics beyond the Standard Model (e.g. from FCNC) or any contribution for which the CP dependence 
of the total Hamiltonian does not satisfy the factorization condition. Numerical results have shown 
effects up to 10 % on the neutrino fluxes in the star. The possible impact on observables like the 
electron-fraction, a key parameter for the r-process, requires further investigation. 
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